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Abstract. - It is suggested that the question of existence of a jamming phase transition in a broad class 
of single-lane cellular- automaton traffic models may be studied using a correspondence to the asymmetric 
chipping model. In models where such correspondence is applicable, jamming phase transition does not 
take place. Rather, the system exhibits a smooth crossover between free-flow and jammed states, as the 
car density is increased. 



Traffic flow and the formation of traffic jams have been extensively studied for many years |T|. A 
very useful quantity which has often been used to characterize traffic flow is the relation between the 
density of cars in the road and the traffic throughput. This relation, termed the fundamental diagram, 
was measured empirically in various situations, and was studied in a large variety of models |2 1 1. 
When the density of cars in the system is low, the traffic flow is expected to grow linearly with the 
density. At high densities traffic jams are formed, lowering down the flow sometimes even to a com- 
plete stop (HIS. This enables one to identify three regimes in the density-flow plane: a free-flow 
regime at low densities; a regime of wide moving jams at high densities; and a synchronized flow 
regime, where jams and free-flow coexist, at intermediate densities The question whether the 
transition from one regime to another is a smooth crossover or is a result of a genuine phase transition 
is still not settled in most traffic models. 

In recent years Probabilistic Cellular Automata (CA) models have been introduced to analyze 
traffic flow ||4l|5l[ll. In such models both time and space are discrete and the physical state of the 
system (e.g. position and velocity of all cars) is updated simultaneously according to some update 
scheme. This provides a rather efficient way for carrying out numerical studies of the fundamental 
diagram. Some traffic CA models were suggested to exhibit phase transitions Il6ll7ll8ll9l . However, the 
existence of such a transition can only be explicitly demonstrated in some limiting cases where certain 
dynamical processes are deterministic. The existence of a jamming phase transitions in more generic 
cases, where all dynamical processes are non-deterministic, is still an open question. Such transitions 
have been suggested to occur in some models on the basis of mean-field methods and numerical 
simulations, which cannot yield a definitive answer to this question. In this Letter we address this 
issue in more detail. 

Within a more general framework, the existence of phase transitions in one-dimensional driven 
systems has been studied rather extensively in recent years and several mechanisms for such transitions 
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have been proposed. Examples include the zero range process llOllll . two species driven models |[T2| 
and the chipping model (2J[l4|. The chipping model introduced by Majundar et. ah incorporates 
dynamical processes which are closely related to those taking place in traffic dynamics, so one would 
hope to obtain useful insights for traffic jams from what is known about the chipping model. 

In this Letter we examine the correspondence between traffic CA models and the chipping model 
more closely. This correspondence suggests that for a large class of traffic models with non-deterministic 
dynamics, a genuine phase transition is not expected. Rather, these models exhibit a smooth crossover 
between the free flow and the jammed phases. In the following we briefly review the main known 
results for the chipping model. We then introduce a simple CA traffic model for which the corre- 
spondence to the chipping model could be made explicit. Other CA traffic models which have been 
introduced and studied in the past are also examined within this approach. 

We start by considering the Chipping Model (CM). The model is defined on a periodic lattice, 
where each site can contain any number of particles. The dynamics is defined through the rates by 
which two nearest neighbor sites containing k and m particles, respectively, exchange particles: 



(fc,m)^(fc + m,0) (/c,m)^(A: + l,m-l) (A:, m) ^(A: - 1, m + 1) . (1) 



The first is a diffusion (or coalescence) process(^) while the last two processes correspond to right and 
left chipping of a particle from one site to the other. It has been shown 1 13 1 that if the chipping process 
is symmetric {ujr = ujl) there is a condensation transition at a critical density, above which one site 
becomes macroscopically occupied. Furthermore, numerical simulations and mean-field studies show 
that the probability P{k) of finding k particles in a site has the asymptotic form P{k) ~ jk^ for 
large k, with r = 5/2. The parameter z < 1 is determined by the average particles density and 
serves as the fugacity. The condensation transition is a result of the fact that r > 2, for which the 
distribution P(k) cannot sustain high densities even atz = 1. This transition is analogous to the Bose- 
Einstein condensation. In contrast, if the chipping is asymmetric there exists no phase transition at 
any density llT4l . In this case numerical studies indicate that the domain size distribution has the same 
form as above, but here r = 2. This distribution remains valid at any density with z approaching 1 at 
high densities, indicating that no condensation transition takes place. 

In the following we argue that the chipping model with an asymmetric chipping process provides a 
framework within which a large class of traffic models can be characterized. Starting from a particular 
traffic model we first identify the domains which characterize the flow. A domain can either be a low 
density segment, termed a gap or a hole in some studies ; a high density segment, termed a jam ; or 
a segment of some other characteristics, defined ad-hoc. A domain of size k is then associated with a 
site of the CM occupied by k particles. One then proceeds by examining the evolution of the domains, 
and identifying their dynamical processes. As will be demonstrated, in many cases these processes are 
closely related to the diffusion and the chipping processes of the asymmetric CM. 

We now consider a particular traffic model for which the correspondence to the CM can be made 
rather explicit. The model, referred as Velocity Dependent Braking (VDB), is a variant of the Nagel- 
Screckenberg model |5|. It is defined on a periodic lattice of size L with M = pL cars. Each car 
is characterized by a velocity Vi{t) = . . . ^max and a position Xi{t). The dynamics is performed in 
parallel by first updating the velocities as 




Xi{t) — 1} with probability 1 — p{vi{t)) 
with probability p{vi {t)) 



(2a) 



and then increasing the position of each car by its speed, 

Xi{t^l) = Xi{t) ^ Vi{t ^ 1) . 



(2b) 



(^)Here only a diffusion to the left is considered. One can also introduce diffusion to the right, without changing the relevant 
results. 
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Fig. 1 - Fundamental diagram for the VDB model. Data for ^ > was obtained from numerical simulations of 
systems of size 500. (a) ^^max = 1 and p = 1/2. (b) ^^max = 2 and p = 0.6 ; here the ^ = free-flow branch is 
given by 2p, while the branch of jammed states in this case was extrapolated from numerical data. 

The braking probability p{v) is defined in terms of the two parameters of the model, p and q, as 



We proceed by first considering the model in the cruise control (CC) limit, q = 0. In this limit 
cars residing in dilute regions of the system move deterministically with the maximal allowed velocity 
^max- It is thus Straightforward to define a free-flow domain as a segment which consists of vacancies 
and deterministically moving cars. For p < pf = l/(vmax + 1) a free-flow steady state exists, where 
all cars move deterministically and the current is simply given by J(p) = Vmax P- At higher densities 
such a state does not exist, local jams are formed and the current is reduced. A phase transition between 
the two regimes thus takes place at some density po < Pf. This particular phase transition is a result 
of the fact that due to the deterministic processes, the free-flow state is an absorbing state which has 
no dynamics. It is thus expected to exist only in the CC limit, and to turn into a smooth crossover for 
q > 0, where no absorbing state exists. 

It is straightforward to analyze the fundamental diagram in the CC limit (q = 0) for the case 
^max = 1- We find Pf = 1/2 and po = {1 — p)/{2 — p). It can be shown that in the jammed state 
(P > Po) the current is J{p) = po{l — p) / {1 — po). For po < p < Pf both free-flow states and jammed 
states coexist in the thermodynamic limit. In this region a free-flow state evolves deterministically and 
jams are never produced. However, starting from a random initial condition a jammed state is formed, 
which slowly evolves towards the free-flow one. However, the time it takes for a system to reach a 
free-flow state increases exponentially with the system size 1 15|. Thus in the thermodynamic limit 
both the free flow and the jammed phases exist as stable steady states of the system. For g > the 
J{p) curve can be calculated numerically, and it shows no singularity, as expected. The fundamental 
diagram of the model for the case Vmax = 1 is given in Fig. [It . 

As in many traffic models, the general features of the model are revealed only at ^max > 1- We 
still expect a genuine phase transition between the free-flow and jammed states in the CC limit. Again 
this transition is expected to turn into a crossover for > (see Fig. [ij)). The question is whether 
there is another transition at p > po which is not associated with the existence of absorbing free-flow 
states, and which may persist beyond the CC limit, namely for > 0. 

In what follows we consider v^ax = 2 in the p > po regime. Here no exact solution is available. 
Instead, we analyze the evolution of free-flow domains, and examine its correspondence to the dynam- 




(3) 



Phase Transitions in Traffic Models 4 




space space 



Fig. 2 - Space-time configurations of the VDB model. Stopped cars are in black, moving cars in gray, vacancies in 
white, (a) The evolution of a single domain is characterized by a biased diffusion and chipping of small domains. 
Every 5th sweep is presented, (b) Chipping and diffusion may lead to coalescence of neighboring domains. 
Configurations are presented every 30 sweeps. 

ical processes of the chipping model. In Fig. |2|we present space-time configurations of the model at 
^max = 2 and p = 0.1. Focusing on a single domain (Fig. |2t) we observe that the domain evolves 
by biased diffusion and chipping processes. In Fig. ^ the evolution of many domains is depicted, 
demonstrating exchange of particles between domains and coalescence of two domains upon contact. 

More explicitly, the dynamics of the domains can be linked to the microscopic processes of the 
traffic model as follows: 

1 . Diffusion & Coalescence - Consider a stationary car, located at the left end of a free-flow do- 
main. This car may accelerate to the maximal velocity Vmax, reach the right boundary and 
eventually stop there. This process corresponds to the diffusion of the whole domain to the left 
(Fig.Et). 

2. Chipping - A stationary car at the left end of a domain may accelerate, and then brake before 
reaching the maximal velocity. This process decreases the size of the domain, and ejects a single 
vacancy to the left (Fig. (S^)). 

Note that chipping occurs only to the left, corresponding to a fully asymmetric CM, with uo^ = {). 

To test this picture we performed Monte-Carlo simulations of the VDB model with q = {) and 
^max = 2, and measured the domain size distribution. For convenience we define the size of a free- 
flow domain, /c, as the number of vacancies it contains(^). We find that the asymptotic form of P{k) 
is consistent with with r = 2 (Fig.|4|), independent of the parameter p. In view of the fact that 



(^)This definition is not unique, as one may include the deterministic cars residing in a domain in the definition of the length, 
without changing the results. One advantage of the definition used here is that the overall length of domains is conserved. 
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Fig. 3 - The Vma^ — 2 dynamics of a domain, depicting (a) diffusion and (b) chipping processes. Stationary cars 
are marked in black and vacancies in white. The velocity of moving cars, marked in gray, is indicated. 
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Fig. 4 - Domain size distribution of VDB model with Vmax = 2, p = 0.4, p = 0.3 and g = 0. Solid line has a 
slope —2. Simulation was performed on a system of size L = 10, 000 and averaged over 15 x 10^ sweeps. 



Fig. 5 - Domain size distribution of TCA model with p = 0.4, a = 0.2, (3 = 0.6, 7 = 0.8 and ^ = 1. Solid line 
has a slope —2. Simulation was performed on a system of size L — 10, 000 and averaged over 60 x 10^ sweeps. 



r = 2, the existence of a macroscopic domain in Fig.|4j as is demonstrated by the peak at large blocks, 
should be interpreted as a finite-size effect, which would disappear in the thermodynamic limit 1 14 1. 

It is interesting to note that unlike the CM, the number of domains in the traffic model is not 
conserved, but is subject to fluctuations. The number of domains is given by M — Mjj where Mjj is 
the number of deterministic cars(^), i.e. cars moving at v = ^max- The fluctuations in Md are expected 
to scale as vT. This was verified by direct numerical studies of systems of size up to 10^. Therefore, 
these fluctuations are not expected to influence the dynamics of the system in the thermodynamic limit. 

We conclude that in the CC limit of the VDB model there is no phase transition for p > po- It is 
interesting to examine the implication of this result on the question of existence of a phase transition 
beyond the CC limit, namely for g > 0. Here the free-flow phase is no longer an absorbing state, and 
thus the transition taking place at po in the CC limit is expected to become a smooth crossover. Since 
no other transition is found at p > po in the CC limit we also expect no such transition for > 0. The 
reason is that in addition to the CM dynamics, a model with q > exhibits other processes, which 
clearly disfavor condensation of macroscopic domains. For example, here a domain may split into 
fractions of comparable size, leading to fast fragmentation of large domains. 

In examining other traffic models, we find that in many cases the dynamical processes character- 
izing the CM may still be used to describe the domain dynamics. However, it may turn out that the 
chipping process involves a detachment of more than a single particle. Nevertheless, as long as the 
number of chipped particles r is bounded by a finite number, or the probability of chipping r parti- 
cles u{r) decays sufficiently fast with r (say exponentially), the main results obtained from the CM 
are expected to be valid. Namely, condensation transition should not take place as long as the chip- 
ping process is asymmetric. To verify this point we studied numerically chipping models for the case 
u{r) ^ exp(— r) and for the case where r is distributed uniformly over a finite range. In both cases 
we found that the domain size distribution decays as z^/k'^ for large k, as expected from the CM. The 

(^) Alternatively, one may associate each domain with a stationary car, say, the one to its left. Within this definition the 
number of domains is given by the number of stationary cars, where domains of size zero are counted as well. Such a domain 
corresponds to two adjacent stationary cars. For the purpose of the present study each definition can be used. 
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VDB model with ^max > 2 corresponds to CM with chipping of more than a single particle. However, 
in the CC limit the number of chipped particles is bound by Vmax(^max — Thus the analysis 

presented above for v^ax = 2 remains valid for any finite ^max- 

To demonstrate the more general applicability of the CM picture to traffic models, we briefly 
consider three other traffic models which have been studied in the past. 

Model for Emergent Traffic Jams 1 16 1 — This is another variant of the Nagel-Screckenberg traffic 
model. The model is closely related to the CC limit of the VDB model defined above, except that the 
velocity update rule (Eq. |2a|l is replaced by 



Unlike the VDB model, here non-deterministic cars slow down rather then brake, and thus with Vmax = 
1 the two models are identical. This model was shown to exhibit self-organized criticality. Based 
on general considerations and numerical simulations 1 16|, it was argued that the size distribution of 
the domains (termed holes in [16]) behaves asymptotically as Indeed, in |16| the evolution of 
domains was described in terms similar to those of the CM. Note that in this case the chipping size 
distribution u{r) decays exponentially with r. 

Traffic Cellular Automata 1 17| — This model belongs to a different class of CA traffic models 
where no velocity variable is attached to a car. Cars move to their nearest neighbor site with a proba- 
bility that depends on the configuration of their neighborhood. In f\l\ the dynamics is defined as 

ot (3 5 
• • oo > • o mo o • o« > o o •• m • om > • o •• o • oo > o o , (5) 

where • denotes a car and o a vacancy. In the symmetric CC case, 7 = J = 1, the model exhibits 
a low-density absorbing state at p < 1/3 and a high-density absorbing state at p > 2/3. It has 
been suggested UTI that for intermediate densities (1/3 < p < 2/3) and in some region of the 
a, /3-plane, the system exhibits a macroscopic jam, suggesting a jamming phase transition at some 
density. The correspondence between this case and the CM is less transparent, and will be addressed 
in a future publication 1 15|. Here we consider the CC limit, 5 = 1, with 7 < 1, and apply the CM 
approach to analyze the jammed phase. In this phase a typical microscopic configuration is given by 
an alternating left-to-right sequence of (a) free-flow regions, composed of vacancies and cars separated 
from their nearest neighbor cars by at least two vacancies, (b) finite mixed region of alternating cars 
and vacancies, and (c) an uninterrupted sequence of cars. In order to apply the approach described 
above, it is convenient to define a domain as a union of adjacent free-flow and mixed regions (a 
and b above). By examining the dynamics of such domains one finds that they indeed exhibit the 
characteristic processes of the CM 1 15 1. We performed Monte-Carlo simulations of this model in the 
CC limit and measured the domain size distribution (Fig.|5l). We find that the distribution is consistent 
with as expected from the CM picture, indicating again that a phase transition does not take place 
in the jammed state. 

It is interesting to note that if instead of (5 = 1 we consider 7 = 1, the role played by cars and va- 
cancies is interchanged. Here one can define a domain as a stretch of cars, within which deterministic 
vacancies are embedded, followed to the left by a mixed region. In this case the distribution of these 
domains (or jams) behaves as and thus no macroscopic jam is expected. 

Traffic Model with Passing — Unlike all models mentioned above, this model is not a cellular 
automaton. The model is defined on a continuous ring, and evolves in continuous time. Each car is 
assigned a-priory a random velocity with which it moves on the ring. When a car encounters a slower 




miIi{Vi{t) + 1, Vmax,^i+l(^) - Xi{t) - 1} 




f \) with probability 1 — p{vi{t)) 

{0, Vi{t + ^) - 1} with probability p(vi(t)) 



(4) 
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car it assumes its velocity, thus creating a jam. With some finite probabiHty the next to leading car 
in a jam can bypass its predecessor and recover its original velocity. The direct correspondence with 
the CM was already noted in iTs"]. Treating the CM within mean-field approximation, the authors 
concluded that the model should exhibit macroscopic jam. However, numerical simulations show that 
the jam size distributions is again k~'^ fTsl . Our approach suggests that in the observed distribution is 
related to the asymmetric nature of the chipping, a feature which cannot be captured in mean-field. 

In summary, it is suggested that in many traffic models the coarse-grained dynamics of domains (of 
either high or low density) in some deterministic limit, may be described by the two basic processes 
of the chipping-model, namely diffusion and chipping. Analysis of several traffic models within this 
approach indicates that as in the asymmetric CM, the traffic models do not exhibit a jamming transition 
beyond perhaps the one related to the existence of an absorbing state. It is concluded that in non- 
deterministic traffic models jamming phase transitions do not take place. Rather, a smooth crossover 
between a free-flow and a jammed state takes place as the car density is increased. The approach 
outlined in this paper could provide a useful tool for analyzing the behavior of traffic models. In 
studying a specific model one first has to establish (using numerical or other methods) that indeed the 
coarse-grained dynamics of the domains does follow the basic processes of the chipping model. Only 
then one can apply the correspondence between the two. It would be of interest to test the applicability 
of this approach to broader classes of traffic models. 

* * * 

We thank M.R. Evans and J.L. Lebowitz for useful discussions. The support of the Israeli Science 
Foundation is gratefully acknowledged. 
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